Effective refractive index tensor for weak-field 
gravity 



Petarpa Boonserm f, Celine Cattoen f, Tristan Faber §, 

Matt Visser ||, and Silke Weinfurtner % 

School of Mathematics, Statistics, and Computer Science, 
Victoria University of Wellington, 
P.O.Box 600, Wellington, New Zealand 

Abstract. 

Gravitational lensing in a weak but otherwise arbitrary gravitational field can be 
described in terms of a 3 x 3 tensor, the "effective refractive index" . If the sources 
generating the gravitational field all have small internal fluxes, stresses, and pressures, 
then this tensor is automatically isotropic and the "effective refractive index" is simply 
a scalar that can be determined in terms of a classic result involving the Newtonian 
gravitational potential. In contrast if anisotropic stresses are ever important then the 
gravitational field acts similarly to an anisotropic crystal. We derive simple formulae 
for the refractive index tensor, and indicate some situations in which this will be 
important. 
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1. Introduction 

Weak-field gravity in Einstein's general relativity is actually more general than 
straightforward Newtonian gravity PQI21- While the approximate validity of Newtonian 
gravity is certainly limited to the weak-field regime, Newtonian gravity makes significant 
additional assumptions as to the smallness of effects that depend on the internal 
stresses, pressures, and energy fluxes in the massive bodies that act as source for the 
gravitational field. While there is no significant doubt that for planets, and most stars, 
the gravitational effects of internal stresses can safely be neglected, the situation for 
neutron stars (where G^M/R « 1/10) is much more uncertain. Furthermore, while 
there is little doubt that the "dark matter" that makes up approximately 90% of most 
spiral galaxies can be treated using weak-field gravity, in the absence of solid physical 
motivation for some particular equation of state we cannot necessarily conclude that 
the gravitational field can be adequately described by Newtonian gravity. 

In view of this we have developed a formalism that makes no assumptions about 
the relative smallness or isotropy of internal stresses (and pressures and fluxes), to see 
how gravitational lensing is affected. In particular, weak Newtonian gravitational lenses 
can be interpreted in terms of an analogy wherein a gravitational field is assigned an 
"effective refractive index" [3J H] , and we extend these ideas to see how this "effective 
refractive index" is affected by the presence of significant internal stress. Most strikingly 
we will see that the "effective refractive index" is in general no longer a scalar, but is 
instead a 3 x 3 tensor — in analogy to the situation in an anisotropic crystal. (The use 
of analogies to relate otherwise distinct phenomena, and to give qualitative insight as to 
what physical effects might be important, has recently attracted significant interest in 
the general relativity community [Sj, but related ideas under the name "electro-optical 
analogy" have an independent history jH].) We organize the paper as follows: 

• First we consider the static case where there are no internal energy fluxes (so in 
particular we neglect the effects of rotation). 

• Second we further specialise this discussion to situations of static spherical 
symmetry - - believed to be a good approximation for galactic halos containing 
"dark matter" . 

• Third we extend the discussion to the more general stationary case, where internal 
fluxes are included. 

• Fourth we indicate how time- dependent situations can in principle be dealt with. 

• Fifth we indicate how the present formalism matches to the usual idea of a far-field 
multipole expansion. 

• Finally we briefly discuss astrophysical situations in which the issues raised in this 
article are likely to become important. 
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2. Static case 

For light propagating in curved space along some curve parameterized by A we have 

9ab dX a (X)dX b (X) = 0. (1) 

Looking at the specific case of a weak field, where g ab = n ab + h ab , the gravitational field 
can be considered as a perturbation h ab around the flat space r\ ab . This leads to: 

g ab dX a (X) dX b (X) = r) ab dX a (X) dX\X) + h ab dX a (X) dX b (X) = . (2) 

Here r] ab dX a (X) dX b (X) is no longer zero. For a light ray propagating in a static weak 
field we get 

dX a dX b . 1 , . fdt\ 2 , c , N di ! diJ n 
^^X^X=^ 1 + h - ) [dx) + ^ + ^dAdA = °- (3) 

Choosing X = t equation (JHJ) simplifies to 

dX a dX b 

9ab = (-1 + K) + (Sij + h i:j ) x i x j = . (4) 

We define a "coordinate speed of light" by calculating the norm | \x l \ | of the "coordinate 
velocity of light " x l using dij, the unperturbed background metric for space. This allows 
us to split the velocity into a speed and a direction 



x l = \\x l \\ k l ; \\k % \\ = 1 = V<% k* ki] (5) 

where k is a unit 3- vector. Putting this into equation (j3J), and noting that the h ab are 
small compared to unity we can usefully Taylor series expand, to obtain 



\x 



' '"' - ' l h u -\h ij % i ti + 0{hl h ) (6) 



<k) = th-tt « 1 + 77 h tt + - h {j U k? + 0{h 2 ab ). (7) 



\ 1 • /,•' /,•' 2 " 2 

for the coordinate speed of light. Note that we have adopted units where the physical 
speed of light, c, measured by physical rulers and physical clocks, is always 1. Then the 
spacetime refractive index for light travelling in the direction k is 

1 1 1 , 1 
,— wl + -L + - 

2 2 

We now define the 3x3 refractive index tensor as: 

Uij = i 1 + 2 ^ + 2 ^' ' ^ 
so that n(k) = k % k^ + 0{h 2 ab ). If we adopt the standard definition 

- h ab = h ab_}_ h ^ (g) 

then a brief computation yields 
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To connect these general formulae to the presence of stress-energy, adopt Einstein-Fock- 
de Donder gauge (that is, harmonic quasi-Cartesian coordinates) pQ 

d a h ab = 0, (11) 

and write the Einstein equations in the exact form pQ 

V 2 h ab = -16tt G N T a f , (12) 

where the "effective" stress-energy contains contributions both from "ordinary" stress- 
energy and the "pseudo-energy" of the gravitational field itself: 

T^b = Tab + tab- (13) 

In all situations we are interested in the pseudo-energy is much smaller than the ordinary 
stress energy, nevertheless explicitly keeping the pseudo-energy as part of the effective 
stress energy is a very useful bookkeeping device. (This is an application of the "principle 
of controlled ignorance" espoused in jT].) Due to the assumption that the spacetime is 
static, the effective stress-energy tensor is 

p cS 



rpeS 

1 ab 



We define 



T eS 
<j J 



(14) 



V 2 $ =4nG N fT, (15) 
V%- = 4vr G N Tf, (16) 

where $ is the ordinary Newton potential, but with the effective mass-energy-density 
p cS as a source. The are new post-Newton gravitational potentials arising from 
the effective internal pressures and stresses of the source matter. Using the Einstein 
equations we find: 

V 2 /^- = - 16tt G N Tf. (17) 
V 2 h tt = - 16tt G n p cS . (18) 

Imposing suitable boundary conditions at spatial infinity and working in terms of the 
density and pressure potentials, we obtain 

htt = -4$; (19) 
hy = - 4%, (20) 

and the equivalent (though more complicated) 

h a = -2($ + 6 kl * kl ); (21) 

hij = - 2(2% + - 5 kl *«]). (22) 

Thus, the refractive index tensor defined in (|8|l()jl takes the particularly simple form 

n l3 = (1 - 2$) <5y - 2%j . (23) 

In view of our use of the effective stress-energy, these statements are now all exact, 
but formal. Because h a b and its derivatives occur in the pseudo-tensor t ab , the RHS of 
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these formulae depend implicitly on h ab - In practical calculations one would start by 
setting h ab — > on the RHS, and then iterate the equations to obtain the desired level 
of accuracy 

For instance, if the effective internal stresses are isotropic Tf^ — > p eS 5{j and 

— > \I/o Siji w e have the simplification 

ha = -2($ + 3^ ); (24) 
h ij = -2^[$-* ]; (25) 
Tly = (1 - 2$ - 2^ ) (26) 

with 

V 2 $ = 4tt Gtv p cff , (27) 
V 2 ^ = 4tt Gjv P eff - (28) 
Thus for isotropic effective stress, the refractive index is a scalar. If furthermore the 
gravitational field is weak and smooth enough that t a b <C T ab this can be turned into an 
approximate statement about physical perfect fluids, where T^- — ► p Sij — the refractive 
index would now be approximately isotropic. (This is the usual situation in most stars: 
the stellar material is to a good approximation a perfect fluid and the gravitational 
pseudo-energy is a small faction of the total mass budget.) 

Finally, if the internal stresses are completely negligible (which is the usual situation 
in most asteroids, but not stars and planets [1.) we can naively set — > to obtain 
the Newtonian limit and recover the well-known standard results [HUE]: 

htt = - 2$; hij = -2$ 8^ riij = (1 - 2$) . (29) 

The novelty in the current analysis lies exactly in the manner in which internal stresses 
in the body generating the gravitational field lead to a "stress potential" which 
then influences both the weak-field metric ()19H22j) and the effective refractive index 
tensor 



3. Static spherically symmetric weak field 

These general considerations can be made more explicit by working in situations of 
spherical symmetry. Consider now a general static spherically symmetric weak field, for 
which the stress energy tensor (with the background Minkowski metric 7] a f, written in 
terms of spherical polar coordinates) takes the form 

" p cfr (r) 

pf{r) 
pf(r)r 2 

pf(r)r 2 sin 2 9 

Here p^ ff and pf ff are the radial and transverse pressures respectively. At the origin we 
must have p° S (0) = p1 s (0). Because of the way the Einstein equations interact with the 
harmonic gauge condition, is exactly conserved in the flat space covariant sense: 

[T^] :b = 0; implying T ab .. b = -t ab .. b (31) 



rpeS 

1 ab 



(30) 
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where the colon denotes the flat-space covariant derivative corresponding to r]^ in 
spherical polar coordinates. Consequently 

2{pf(r) -pf(r)} 



d rP f(r) + 

T 

which we can use to eliminate pt{r) as 



0. 



pf(r)= P f(r) + -rd r pf(r). 



(32) 



(33) 



We can similarly write 



h, 



ab 



H (r) 
H r (r) 
H t {r)r 2 









(34) 



H t (r)r 2 sin 2 6 

where at the origin symmetry demands H r (0) = H t (0), and the harmonic gauge 
condition (fTTj) . now in the flat-space spherical polar sense h ab -,h = 0, implies 



H t (r) = H r (r) + -rd r H r (r). 



(35) 



After a brief computation 



V 2 h 



ab 



h :c 
,L ab:c 



Z (r) 

Z r (r) 










Z t (r) r 2 

ZAr)r 2 sin 2 6 



(36) 



where we are again using the covariant derivative for flat-space spherical polar 
coordinates, and we have defined 



Zo = ^d r (r 2 d r H ) ; 
Z r = -^d r (r 2 d r H r 



Z t = \d r (r 2 d r H t ) + 2 



HJr) - HAr) 



HJr)-H t (r) 



Z r + — t d r Z r . 



(37) 
(38) 
(39) 



The slightly unusual terms proportional to the difference H r — H t arise due to the fact 
that we are now using spherical polar coordinates (not strictly harmonic coordinates), 
and because has tensor indices in its own right. The Einstein equations now yield 

-16tt G N p eS - 



Z 



-167r G N p e r 



eff . 



plus the redundant equation 

Z t = -16tt G N pf. 
These have the formal solutions 



H = -16?r G N J (J p eff r 2 dr|r- 2 dr = -4$; 



(40) 
(41) 

(42) 
(43) 
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H r = -16tt G N J y pf r 4 dr| r" 4 dr; (44) 

where we have again used equation (JTHJ), and so 

H r -H t = 8 -^{J P fr*dr}. (45) 

We can further rearrange this result slightly by defining 

= H r + 2H t ; H A = H r - H t , (46) 

so that 

H, = ^±^; H, = (47) 

For a null curve located at radius r, and making an angle \ with respect to the f 
direction, equation (fTUJ) implies 



so that 



or 



n{r, X ) = l + ^H + - (H r cos 2 X + H t sin 2 X ) , (48) 
n(r, X ) = 1 + Ih + + ^ A {3 cos 2 X ~ l} , (49) 



n(r, X ) = 1 - 2$ + ltf s + l -H A {3 cos 2 X - l} ■ (50) 

This particularises equation ()23|) to the case of spherical symmetry, with now being 
written in terms of and and with explicit integral formulae now being available 
for these quantities. 

Let's compare this with a situation where we know the exact result: If the object 
that acts as source for the gravitational field has a definite surface, with vacuum outside 
that surface, then application of the Birkhoff theorem tells us that the spacetime 
geometry will be Schwarzschild outside that surface. But because of our subsidiary 
assumption, that we are working in the spherical polar version of Einstein-Fock- 
de Donder gauge, there is no remaining freedom in our coordinate system and we must 
obtain the weak-field limit of Schwarzschild geometry in harmonic coordinates — for 
which the effective refractive index is asymptotically isotropic. In harmonic coordinates 
the Schwarzschild solution has 

h rr=- - + 0{r-')- (51) 



2m m 2 



/■-' r 2 sin 2 9 r r 2 



+ -?\ (52) 



H A = - 3 -^ + 0(r~ s ). (53) 



Hence anisotropies in the effective refractive index are (at least in the case of spherical 
symmetry) constrained to rapidly decay outside the source that generates them. We 
now show that this is compatible with equation (|43j) above. Outside the body T ab = so 
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pf 1 — ► ^rr depends on the pseudo-stress only. But t rr ~ (dh) 2 ~ (m/r 2 ) 2 = m 2 /r 4 . This 
then implies if a ~ m 2 /r 2 , in agreement with the exact result. The precise distribution 
of effective refractive index anisotropies inside the source object can only be determined 
by solving the Einstein equations, which requires some specific model for the interior 
distribution of anisotropic stresses. This is a topic which we hope to explore more fully 
in the future. 



4. Stationary case 

Let us now consider the stationary case, where 

htt h t j 



h 



ab 



(54) 



and due to symmetry htj = hjt- The condition for a photon trajectory then becomes: 
dX a dX b . 1 , . f dt\ 2 oL dt dx* , s dx*dxj n 



Choosing the parameter A = t, this becomes a quadratic equation in x l with the 
coordinate speed of light now being given by 

ii^ii = - — ^-frjr (- htj y + ^(h t] yy-{-i + h tt ){5 tJ + h tJ )yy) . (56) 

[pij + hij) k l k3 V / 
Simplifying, Taylor expanding, and inverting gives the refractive index 

n (k) = i + ~ (h tt + 2^ y + h^ y y) + o(h 2 ), (57) 



2 

which is a very straightforward extension of the static case. But because of the linear 
term in k, it is not possible to bring this completely into 3-tensor form — there is 
additional structure and we must write 

n(k) = Hij y y + h tj y + o{h 2 ), (58) 

where has exactly the same form as in the static case [see equation (jHJ)] and 
the new h t j A/ 7 term can be interpreted as being due to motion of the "effective 
medium" with respect to the quasi- Cartesian coordinate system (t, x l ). This can be 
justified by performing a coordinate transformation into the local rest frame of the 
"effective medium", which is moving with "velocity" —htj, and exhibits a refractive 
index tensor n^. See Appendix A for details. Since the "medium" is generally moving 



inhomogeneously the local rest frame is not best for performing explicit calculations 
- for practical calculations it is preferable to fix the coordinate system once and for 
all and to work with the n(k) of equation ()58|) above. (In particular, the differential 
equations relating the weak field h ab to the distribution of stress energy are defined in the 
original coordinate system, which was chosen to satisfy the Einstein-Fock-de Donder 
gauge condition, and these differential equations do not adapt nicely to the comoving 
point of view.) 
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We can furthermore define additional "flux potentials" 
V 2 IL 



4tt G N T t f 



(59) 

that couple to the momentum flux. The corresponding weak-field Einstein equations 
are very simple, since f] ab vanishes for the off-diagonal elements: 

V 2 h t3 = V% = -167T G N T t f. (60) 

Imposing appropriate boundary conditions 

hu = -4 IL 



(61) 
(62) 



and the refractive index for the stationary case is 

n(k) = U k j -AUj k? + C($ 2 , ^ 2 , n 2 ) 

where riy($, \1/) is as in equation (j2Sl). 

An alternative representation is to note that if we define k a = (1; k l ) so that k a is 
a null vector with respect to the background metric, then 



n(k) = l + -h ab k a k b + 0(h 2 ) 



1 



1 — 



h ab k a k b + 0(h 2 



If we further define 
t lien 

n{k) - 



1 - 2 $ ab k a k b + o($ 



> 2 a6 ) 



(63) 



(64) 



(65) 



5. Time dependent situations 

In the presence of time- dependent sources the only modification is that the Laplace 
equations for the potentials should be replaced by wave equations V 2 — ► — <9 2 + V 2 , 
which has the effect of replacing the usual 1/r potentials by the appropriate Lienard- 
Wiechert potential jjj. The net result is that 

hab{x,t) = 16ttG n / d y + — — J -, (66) 

J \x — y\ 

where t is the retarded time 

t = t — \x — y\. (67) 

But if the so-called null energy condition [NEC] is satisfied [7j, then in particular 
T ab k a k b > 0, implying h ab k a k b > 0{h 2 ab ), so that n{k) > 1 + 0{h 2 ab ). That is, the NEC 
implies the effective refractive index is [to order 0(h 2 b )] greater than unity, thereby 
guaranteeing that in this approximation the coordinate speed of light is always less 
than 1. This connects the discussion back to the perturbative version of "superluminal 
censorship" discussed in jj]. 
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6. The far field 

For isolated bodies it is possible to expand the far field in terms of multipole moments. 
In particular, if we go to the rest frame of the body, and if the body has time- independent 
internal structure, then it is a standard result that [l] 
2G N M 



lab 



5 ab + 0(l/r 2 ). (6* 



r 

Thus for a single isolated body 

n (k) = l + -h ab k a k b + 0(h 2 ) = 1 + + 0(l/r 2 ), (69) 

2 r 

where k a = (1; & l ) is a null vector with respect to the background metric. But as long as 
the total mass is positive M > we have n(k) > 1. That is, the positive mass theorem 
implies that in the far-field regime the effective refractive index is always greater than 
unity, thereby guaranteeing that the coordinate speed of light is always less than 1. This 
is a rather different perturbative version of "superluminal censorship" , more akin to the 
ideas discussed in jHj. Furthermore note that for an isolated body at rest the far- field 
refractive index is automatically isotropic — this is consistent with our calculation for 
situations possessing spherical symmetry, where we found that anisotropic propagation 
of null geodesies was confined to the region inside the body. 

If there are many compact objects making up the gravitational lens, one should 
simply sum the 1/r potentials for each object. If we are considering a body that is 
at rest but rotating, then the overall momentum P l is zero and the dominant term in 
the potential IT comes from the dipole term in the multipole expansion — this term 
is well-known to be related to the total angular momentum of the rotating body and a 
standard textbook result yields [T] 

u _ G N (Lxv) i + 0{i/r3) (7Q) 

If we consistently retain all 0(l/r 2 ) terms in a post-Newtonian analysis then the far 
field for a stationary rotating body is (see, for example, equation (19.13) of pQ) 



h, 



ab 



" 2G N M/r - 2G 2 N M 2 /r 2 


-2G N (L x i) j/r 2 


-2G w (Lxf),/r 2 


(2G N M/r + 3G 2 N M 2 /2r 2 )5ij _ 



+ 0(l/r 3 ) (71) 



which implies 

. = l + gg^M _ _ 2G, (L x r) ■ k + 



So angular momentum certainly does leave an imprint on the far-field refractive index, 
but in a sub-dominant term. 

Of course if the bodies in question are not isolated (one might for instance be 
considering a null geodesic passing through the bulk of the body), then the multipole 
expansion is not useful — and one should resort to the use of the general equation (jH2j) 
and the linearized field equations for $, 11*, and (with V 2 — > — d 2 + V 2 for time 
dependent situations). 
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7. Discussion 

Our formulation of weak field [but not necessarily Newtonian] gravitational lensing is 
particularly simple and gives a nice interpretation of the "effective refractive index" 
of the gravitational field directly in terms of Newtonian-like potentials coupled to the 
stress-energy tensor. Doing so results in an "effective refractive index" more general 
than the standard Newtonian result - - anisotropic stresses are seen to lead to an 
anisotropic "refractive index" and energy fluxes and angular momentum are seen to 
lead to a "moving medium" effect. 

First we should raise (and settle) the issue of gauge invariance [coordinate 
independence] - our results were obtained by using the simplifying properties of 
Einstein-Fock-de Donder gauge, and so clearly there is a sense in which the existence 
and value of the refractive index tensor depends critically on the use of specific 
coordinates. However it must be emphasised that once one has used the refractive 
index to calculate quantities such as magnification factors and/or the angle measured 
between two images on the sky, these angles and magnifications have a coordinate- 
invariant physical meaning independent of whatever coordinate system was used to 
carry out the calculation. So while the analogy that leads to the refractive index tensor 
is not itself coordinate invariant, the physical observables that result at the end of any 
specific calculation are coordinate invariants. 

We now ask under what circumstances these effects might physically be important? 
Are there physical situations in which the gravitational field is appreciable (but still 
a weak-field), but with « T tt l Perhaps the best known situation of this type 
arises in the core of a neutron star where GM/R « 1/10 so the gravitational field 
is reasonably weak, and where matter is approximately described by a "stiff" equation 
of state, p = pc 2 . Furthermore, a radiation fluid (or neutrino fluid) satisfies p = pc 2 /3, 
so if radiation fluid ever becomes a significant fragment of the total mass budget of any 
clumped system, the ideas of this article would become important. The contribution 
of angular momentum to microlensing, which in our formalism arises at the dipole 
level in the Ilj potential, has been considered in several articles [5] using somewhat 
different formalism. A particularly interesting effect is the possibility of anisotropy in the 
refractive index — this occurs once one moves away from the perfect fluid approximation 
for astrophysical bodies, and speculation concerning significant crustal stresses in 
neutron stars (and their possible effects on compactness bounds) is common ^Hj- Finally 
we point out that coherent field configurations, be they electromagnetic fields or scalar 
fields, generally induce anisotropic stresses comparable in magnitude to the energy 
density. In particular any attempt at modelling the "dark matter" in galactic halos 
with classical fields will lead to anisotropic stresses that, while weak, are comparable in 
magnitude to the energy density [TTj. 

In summary, so long as there is continuing uncertainty over the total mass budget 
and relevant equations of state for the various compact objects occurring in our universe 
we feel that it is prudent to retain the generality of the analysis in the present article. 
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To see how a stationary (non-static) weak gravitational field can be viewed as a moving 
"effective medium", it is useful to temporarily work with strong gravitational fields. 
Pick a particular point in spacetime (to,x l ). and consider the change of coordinates 

*->* = *; x i -> x l = x i -v\tQ,x Q ) [t -t ] (A.l) 
with v l (t ,x ) = g lk (t ,x ) g t k(t ,x ). Then 

dt = di; dx l = dx* + g ik (t , x ) g tk (to, x ) dt (A.2) 
and in the immediate vicinity of (to,x l ) we have 

g ab dX a dX b = g tt dt 2 + 2g ti dt dx l + 9ij dx' dx j (A.3) 

= [gu - (gu g ij g^W + gij dx* dx j (A A) 

which clearly has the effect of (locally) banishing the mixed time-space parts of the 
metric at the cost that 

gu = gu - (gu g 13 gjt); gu = 0; g tj = g^. (A.5) 

But because we wish to apply this in the weak field approximation g ab = i] ab + h ab with 
h ab 1, this simplifies tremendously Since g tl = h tl + 0(h 2 ), whereas g ti = h ti exactly, 
we see 

h tt = h tt + 0(h 2 ); h a = 0; hj = K- ( A - 6 ) 
That is, in weak field going to the moving coordinates defined by 

t->t = f; x i ^x l = x l + 5 ik h tk (t ,x ) \t -t ] + 0(h 2 ); (A.7) 
dt = dt; dx* = dx l + 5 ik h tk (t , x ) dt + 0(h 2 ); (A.8) 

allows us (locally) to interpret — h tk as the "velocity" of the medium, while and h a , 
since they change by at worst 0(h 2 ), lead to the same refractive index tensor as was 
encountered in the static case. 
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